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Abstract 

We recently proposed a new approach to the Casimir effect based on classical ray optics (the 
"optical approximation"). In this paper we show how to use it to calculate the local observables 



We work three examples in detail: parallel plates, the Casimir pendulum and a sphere opposite 



■ of the field theory. In particular we study the energy-momentum tensor and the Casimir pressure. 
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f-*) I a plate. We also show how to calculate thermal corrections, proving that the high temperature 
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■ 'classical limit' is indeed valid for any smooth geometry. 
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I. INTRODUCTION 

The Casimir effect is a manifestation of the quantum fluctuations of a 

quantum field at a macroscopic level. Experiments on Casimir forces are precise tests of one 
of the less intuitive predictions of field theory. For a theoretician, predicting the outcome of 
these experiments is a worthy challenge. Hence it seems somewhat astonishing that an exact 
solution exists only for infinite, parallel plates case 1]. Other formal solutions for geometries 
not made of distinct rigid bodies free to move (like the wedge, the interior of a sphere 
or of a rectangular box 3]) are irrelevant for an experimental setup. Moreover in many 
such solutions divergences have been discarded in a way that leaves the result unrelated to 

n 

practical materials and configurations [6]. Interesting theoretical developments include the 
method developed in [3| where the solution for infinite periodic geometries is obtained as a 
series expansion in the corrugation, and the numerical Montecarlo analysis in ^. 

Amongst the various effects that an experimentalist must take in account to interpret 
the data (e.g. finite conductivity, temperature and roughness corrections) probably the most 
challenging, interesting and full of connections with other branches of physics and mathe- 
matics is the dependence of the force on the geometry of the bodies. Calculating the Casimir 
force for perfectly refiecting bodies in the end reduces to finding the density of states (DOS) 
of the Scrhodinger Hamiltonian for the equivalent billiard problem including the oscillatory 
ripple on the averaged DOS. This is an incredibly difficult problem in spectral theory that 
still challenges mathematicians and physicists today ^ and in essence is not solved beyond 
the semiclassical approximation. 



In this context we have introduced in Refs. [10|, ^| a method based on classical optics 
which has several virtues: accuracy, uniform validity when a symmetry is born, straightfor- 
ward extension to higher spin fields, to non-zero temperatures, to include finite refiectivity 
and, the main topic of this paper, it can provide an approximation to local observables. 

This paper is structured as follows: in Section |n] we show how to cast the energy mo- 
mentum tensor into a sum over optical paths contributions and how to regulate and analyze 
the divergences, ubiquitous in Casimir energy calculations. Section III is dedicated to the 
analysis of the three examples already studied in with pedagogical intent. We study 
parallel plates, the Casimir torsion pendulum and a sphere opposite a plate. In Section IV 
we show how to calculate the same local observables and the free energy for a thermal state 
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and we prove (within the hmits of our approximation) the 'classical limit' theorem |4l Il2j|. 
which states that at high T, Casimir forces become independent of h and proportional to 
T. As far as we know this is the first time this assertion can be generalized to geometries 
other than parallel plates. We also study the example of parallel plates (finding the known 
results) and of a sphere opposite a plate at non-zero temperature. We find evidence, again 
within the framework of the optical approximation, that the low T behavior of the Casimir 
force is a difficult problem, qualitatively different from the T = and high temperature 
cases. 



II. LOCAL OBSERVABLES 



Local properties of the quantum vacuum induced by the presence of boundaries are of 
broad interest in quantum field theory Q]. For example gravity couples locally to the energy- 
momentum tensor. Vacuum polarization induces local charge densities near boundaries, 
provided the symmetries of the theory allow it. Also, local densities are free from some 
of the cutoff dependencies that plague many other Casimir effects. Any local observable 
that can be expressed in terms of the Greens function can be estimated using the optical 
approach. In this section we study the energy, momentum and stress densities for a scalar 
field. 

Some local observables are not unambiguously defined Q|. For example the charge 
density (in a theory with a conserved charge) is unambiguously defined while the energy 
density, in general, is not (while its integral over the volume, the total energy, is). In this 
paper we use the Noether definition of the energy-momentum tensor, similar results would 
be obtained with other interesting definitions. 



A. Energy-momentum tensor 

We study the Noether energy-momentum tensor of a free, real scalar field in a domain V 
with Dirichlet boundary conditions (BC) on S = &D made of (in general disconnected) sur- 
faces. Other BC (Neumann, Robin) can be discussed but for simplicity we restrict ourselves 
to Dirichlet BC here. 
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The lagrangian is (we use h — c— 1) 

C^]^d^4>d^<t^-\m^\ (II.l) 

where Greek letters are used for 4-dimensional indices while the vector notation will be used 
for spatial vectors. 

The Noether energy-momentum tensor for this real scalar field is 

dC 

T^. = - g,r\ {da<f> - (II.3) 

from which we identify the energy density Tqq, the momentum density Toj, and the stress 
tensor Tij. The definition of these quadratic operators involves divergences that wc will regu- 
late by point splitting. We hence replace quadratic operators like 0(a;)^ by lim^'^x (p{x')(p{x). 
The energy density operator, for example, is 



Too{x,t) = lim 



^do(p{x', t)do(i>{x, t) + ^V' • V(j){x', t)4>{x, t) + ]^m^{x\ t)(j>{x, t) 
^do(i){x',t)do(i){x,t) - ^(j){x',t)V^(j){x,t) + 



= lim 

+\m^{x\ t)<f>{x, t) + ^ (V' + V) • t)V<f>{x, t)] (II.4) 
The field satisfies the free wave equation in T> 

d"^^ + w?4) = (II.5) 
and hence it can be decomposed into normal modes 

'^(^' ^) = E 7ff = + ^;(:^^)e'^^*a]) , (11-6) 
where if^j and Ej are the eigenfunctions and eigenvalues of the problem 

{-V'^ + m'^)^j = E]^j forxeP; i^jix) ^ for x e <S. (11.7) 

We also use the definition E{k) — -\/P~+m^, and Ej — y^fc|~+m^ so that the eigenvalue 
equation reads 

- V% = (11.8) 
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and because of the positivity of the operator — V^, the spectrum {Ej} is contained in the 
half-hne {E > m}. 

We now introduce the propagator G{x',x,k), defined as in Ref. to be the Green's 
function of the problem ()II.7|) or ()II.8|) : 



{-W'^ -k^)G{x',x,k) = 6{x'-x) 

G{x', x) = for x' OT xe S, (11.9) 

which can be written using the spectral decomposition as 

In Ref. j3] we have developed an approximation for the propagator G{x', x, k) in terms of 
optical paths (closed, in the limit x' —>■ x). The derivation can be found in Ref. ll|, the 
general result valid for spatial dimensions being 

Gopt{x',x,k) = r^lv/9.4-^^^r/■?.^^ (ir^r) ^ k^/"^ ^H^^^_^{k£r) , 

r 

= J2'^r{x',X,k), (11.11) 
r 

where if is a Hankel function, r labels the paths from x to x', Ur is the number of reflections 

of the path its length and Ar{x',x) is the enlargement factor familiar from 
classical optics, 

Ar{x',x) = f^. (11.12) 

Ar{x', x) is the ratio between the angular opening of a pencil of rays at the point x and the 
area spanned at the final point x' following the path r. For = 3 we have 

GJx',X,k) = (-l-^n Ar^^(^'^^) ^ikM.',.)_ .jj_^3) 

An 

With this explicit form for the propagator G, we now have to rewrite the elements of the 
quadratic operator T^j, as functions of G and its derivatives. It is useful to pass from the 
point-splitting to a frequency cutoff by inserting the latter in the normal modes decompo- 
sition pi.6|l as 

-fcj/A 



e 



dke-^'^2k5{k'' -k]). (11.14) 
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The limit x' ^ x can then be exchanged with the dk integral and we get for the energy 
density, 

/OO 1 /.OO 7 

dke-'/''-E{k)p{x,k) + rfA;e-^/A__v./(a;,A;). (11.15) 

The density p and the vector j are defined as 

2k 

p{x,k) = — lmG{x,x,k) (11.16) 

1 - , 1 - 

j{x,k) = lim — ImVG(x', X, A;) = — ImV G{x,x,k). (H-l^) 

x'—^x TC 2tI 

£ is obtained by integrating Too over the whole volume between the bodies: 

£= £x dke-^'^-E{k)p{x,k)+ / dke-^'^—-— / dS-j{x,k). (11.18) 
Jv Jo 2 Jo 2E[k) jg 

We have turned the integral over the divergence of j into a surface integral using Gauss's 
theorem. In the case of Dirichlet or Neumann boundary conditions, since dS oc n we have 
(here dft = n ■ V and jn = n ■ j) 

jn{x, k) = — Im dfiG{x, x,k) = 0, x E S (11.19) 

71 

and the surface integral term disappears. It should be noted that the vanishing of the j 
contribution to the total energy relies on the continuity of the propagator for x', x E V. 
In some approximations, including the optical one, this continuity is lost. Hence spurious 
surface terms arise on the boundary of certain domains V C V. This region is what in wave 
optics is called the 'penumbra' region. Diffractive contributions are also not negligible in 
this region and they cancel the discontinuities in G, hence eliminating the surface terms. ^ 
The surface terms in the energy are hence of the same order of the diffractive contributions 
which define the error in our approximation. 

The divergence V- j could also be eliminated from Too by changing the energy-momentum 
tensor according to 

f^u = T^u + (11.20) 

with 

Ipat^u = ^(P {gfiuda - Qaudfj) (j). (11.21) 



^ As an example see Kirchoff's treatment of the diffraction from a hole in Ref. [l^. 



The total energy £ and momentum are not affected by this redefinition however the new 
tensor T^^, is not symmetric. 

It can be seen that the stress tensor T^j is normal on the surface S (for both Dirichlet 
and Neumann BC) so locally the force on the surface is given by the pressure alone 

dp 

_ = nP=(0|T^,^|0). (11.22) 
The operator Tfiji regulated by point splitting is 



Tft,fi{x,t) = lim 



— lim 

x'—*x 

- ^ (v' + v) <f)'W(p (11.23) 

where 0' is shorthand for (f){x', t). The second term in brackets is zero when averaged over an 
eigenstate of the number operator by virtue of the equations of motion. For Dirichlet 

BC the term 0V^0 = on the boundaries, so we have (V = ndfi + V*) 

(0|T^,„HO) = ]nn ^ {dLdn - • V, + ^,{^')^Ax). (11.24) 
j ■' 

Since also Vtipj{x) = on the boundaries this expression simplifies to 

P{x) = hm J2 ^m^A^')^A^). (11.25) 

This expression can be rewritten, in terms of the propagator G, regulated by a frequency 
cutoff as we did for Tqo, 

r°° k 

P{x) = lim d'fidfi dke-''/^———lmG{x',x,k). (11.26) 
x'-*x Jq zTrEyk) 

In this regulated expression we can exchange the derivatives, limit and integral safely. Below 
we discuss what the divergences are when A ^ oo and how to interpret and dispose them. 

All the above expressions are exact. Once the propagator G is known, we can calculate 
the energy-momentum tensor components from them. However as discussed above in the 
interesting cases it is difficult to find an exact expression for G and some approximations 
must be used. 
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For smooth impenetrable bodies we use the optical approximation to the propagator 



developed in Ref. 



IG 



111 ] and recalled in eq. This gives G as a series of optical paths 



and hence the pressure P as a sum of optical paths contributions 

P ^ J^Pr (11.27) 

r 

Pr = i-ir lim d'^drt / ^^e-^/^-— -^-p^ sin (H,(x', x)) , (11.28) 

An important feature of the optical approximation is that all divergences are isolated in 
the low reflection terms whose classical path length can vanish as S. In practice 

only the zeroth and first reflection are potentially divergent. Before performing the integral 
in k and taking A ^ oo then we have to put aside the divergent zero and one reflection 
terms Pq and Pi for a moment (in the next section we will show how their contributions are 
to be interpreted). 

For the remaining families of paths (that we will denote as r G 7^) the integral over k 
can be done and the limit A ^ oo taken safely. The result is finite and reads 

We can further simplify this expression. For simplicity let us call z the normal direction. 
Notice that for any sufficiently smooth function f{z', z) vanishing for either z' or z on the 
surface z = Q 

d..dj{z\z) =ld'j{z,z) . (11.30) 

z'=z=0 2 z=0 

The proof is trivial: consider that the lowest order term in the expansion of f{z', z) near 
z',z = is Qc z'z. The propagator G{x', x, k) satisfies all these properties and hence we can 
use this result to get rid of the limit x' x and assume x' = x from the beginning. We can 
therefore rewrite Eq. pi.29|) as, 



Equation (nT3T]) is one of the main results of this paper. In Ref. we reduced the 
computation of Casimir energy to a volume integral. The force is then found by taking the 
derivative with respect to the distance between the bodies. Calculating the pressure instead 
gives the force by means of just a double integral of a local function. The problem is then 
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Ay^{x, x) 



computationally lighter and sometimes (as we will see in the examples) can even lead to 
analytic results. 

Essentially the problem has been reduced to finding the lengths and enlargement factors 
associated with the optical paths for points dose to the boundary. In the case of the pressure 
(eq. pi.29|) or ()11.31|) ) it is necessary to know their derivatives in the direction transverse to 
the surfaces. We will see that this problem can be easily tackled numerically when it cannot 
be solved analytically. 



B. Regulate and eliminate divergences 



As in the energy calculations the only divergences occurring in the pressure come 
from by paths whose lengths i <^ 1/^5 where A is the plasma frequency of the material. 
There are only two such families of paths: the zero and one refiection paths. In this section 
we show that these divergent contributions are independent of the distances between the 
bodies. This fact is easily understood: in order for a path to have arbitrarily small length 
all of its points must be on the same body. So in order to study these terms we need only 
consider a single, isolated body (and a massless field). We are also careful in maintaining 
the double derivative d^, ^ since we are calculating the terms Pq and Pi separately. 

For r = 0, the zero refiection term, introducing an exponential cutoff A on the material 
refiection coefficient we obtain 

n = re-^'-'dk-^ J.na,a. f^ij^g^l = ^. (11,32) 



jQ 2'KE{k) x'^x&s \ Aii\z' - z\ ) Att'^ 

The same calculation for the r = 1 or one refiection term gives: 

= r^-"'dk^— hm d,d. (-'^^^^\ = ^. (11.33) 
Jo 27rE{k) x'^xes "\ A-k\z' + z\ J Air^ ^ ^ 

Notice that these two terms are equal, so we could have substituted d^'^z — ^ for their 
sum, after having properly regulated the divergence. 

This positive, cut-off dependent pressure, P^ = Pq + Pi, must be dynamically balanced 
locally by a pressure generated by the material, lest it collapse. Moreover the total force 
obtained by integrating this quantity over the (closed) surface S of the whole body gives 
zero. However, if the space around the body were inhomogeneous, as in the presence of a 
gravitational field, a finite term survives the surface integration, giving rise to a "vacuum 
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Archimedes effect" in which the pressure on one side is, due to gravitational effects, larger 
than on the other side, so the body feels a net force. We have analyzed this effect in detail 
in Ref. and called it "Casimir buoyancy". 

Finally note that another important element of this class of quadratic operators is the 
Feynman propagator. In studying a field theory in a cavity or in between impenetrable 
bodies (for example hadrons as bags, photons in cavities or Bose-Einstein condensates in 
traps), we can consider expanding the Feynman propagator in a series of classical optical 
paths reflecting off the boundaries. The first term, related to the direct path is the familiar 
free propagator, the others give the finite volume corrections. 



III. EXAMPLES 



In this section we calculate the Casimir force from the pressure, using the formalism 
developed in the previous section, for three examples that were already addressed in Ref. 
[lit using the energy method. 



A. Parallel Plates 



The parallel plates calculation is a classic example, whose result is well known and con- 
stitutes the basis the widely used proximity force approximation (PFA) 3|. We use this 
standard example to establish the rank among contributions to the total pressure and show 
the similarity and differences with the energy method 

We calculate the force acting on the lower plate, denoted by d or down, by calculating 
the pressure on its surface. We discard the zero and Id (one reflection on the lower plate 
itself) reflection terms. The first term to be considered is the path that bounces once on the 
upper plate {u or up) lu. For parallel plates A = and we have 

pi^) = y] i-^Tdi-—^. — -. (III.l) 

The length £^(a;, x) for the paths that bounce an even number of times is a constant in 
z and hence the derivatives vanish: they do not contribute to the pressure. This seemingly 
innocuous observation simplifies the calculations considerably and it is a test for any other 
geometry which reduces to parallel plates in some limit: in this limit the even reflections 
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FIG. 1: Odd reflection paths that contribute to the Casimir force between the two plates in the 
pressure calculations with the optical approximation. The points x' and x will eventually be taken 
coincident and lying on the lower plate. 

contributions must vanish. Generically their contributions are small. This parallels the role 
of the odd reflection paths in the energy method 11 1. 

Figure Q shows the odd reflection paths labelled with our conventions. For the path \u 
we have 



327r2a4 

The next path to be considered is the path that bounces 3 times, first on (i, then on u 
and again on ci, dud = 3u (3 stands for 3 reflections and u for the plate where the middle 
reflection occurs) which gives a contribution 

1 3 



P^uix) = \im—d^ — -7- 



(III.3) 



327r2a4 

The two contributions Eq. pil.2j) and Eq. pil.3j) are equal. The reason is easily uncovered. 
One can recover Eq. (fnT3|l from Eq. (fnT2|) sending z —z but for the purpose of taking 
the second derivative at z = this is irrelevant. In the same fashion = P-^^i P^u = Pju 
etc. and hence we find 

,4 



P(x) = -2- 



+ ... 



(III.4) 



327r2a4 327r2(2a)4 327r2(3a)4 167r2a4 90 

which is the well-known result. Notice also that the rate of convergence is the same as in 



the calculation making use of the Casimir energy in Ref. [ll| (?7,-th term contributes 
of the first term, in this case lu + 3^). These observations that allow us to determine the 
rank of the contributions are fundamental, and they apply as well to the other examples in 
this section. 
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B. The Casimir Torsion Pendulum 

In this section we study a geometry already considered in Ref. a plate inclined at 
an angle 9 above another infinite plate. We have called this configuration a 'Casimir torsion 
pendulum' because the Casimir force will generate a torque which can be experimentally 
measured. The configuration is analogous to the parallel plates case but the upper plate 
must be considered tilted at an angle 9 from the horizontal. The length of the upper plate 
must be taken finite, we denote it by w, while the length of the lower plate can be infinite 
which we choose for simplicity. There is only one substantial difference with the parallel 
plates case: the even reflection paths do contribute in the pendulum, since their length varies 
as we move the flnal points x', x. 

We calculate the force exerted on the lower, inflnite plate for simplicity. We then obtain 
the energy by integrating over the distance along the normal to the lower plate and from 
this we can calculate the torque as 

r^-|. (....5) 

The lower plate is taken inflnite, the upper plate width is w, and the distance between 
the height at the midpoint of the upper plate is a. We will choose as the origin of the 
coordinates one point on the intersection line between the lower plate and the line obtained 
by prolonging the upper plate. This deflnes a flctitious wedge of opening angle 9. We call 
X the horizontal and z the vertical coordinate, the third direction, along which one has 
translational symmetry, being y. 

Since the surfaces are locally flat we have A = as in the case of the parallel plates, 
and again the odd reflections are exactly as in the case of the parallel plates. However now 
the even reflections contribute (the notation is the same as in the parallel plates case, in the 
even reflections 2u means the flrst reflection is on the upper plate etc.): 

P = Plu+3u + P2u+2d + -Psd+Sd 

+ ... (III.6) 

where we have grouped the terms with the symbolic notation Pa+b = Pa + Pb when Pa = Pb- 
It is useful to recapitulate what we have learned about the rank of these contributions: 
Piu+3u dominates, P3d+5d is smaller by ~ 1/16, P5u+7u is smaller by ~ 1/81, etc. The even 
reflections are generically much smaller than the odd reflections, and vanish as 6* 0. 



12 



The first term in ()III.6|) is 

1 2 1 

Plu+3u = —2 2 ^2^27 V' (J^^-'^) 

with £i = 2{xsm6 — zj cos 6), and an overall factor of 2 takes into account the identity 
Piu = P-su- Taking the derivative and then setting z = we find 

3 1 

Plu+3u — ~T7. 2^ '■ 47) V7i'> (III. 8) 

Ibvr^x^ sm 6'cos"'^ 

and integrating from Xm = {o-/ sai9 — w/2)/ cos 6 to xm = {a-/ sinO + w/2)/ cos 9 we find the 
force per unit length in the y direction 



_ cos^ / 1 1 



fill 9) 

32n^sm^9 \{a/sme {a/ sind + w/2y / ' ^ 
Since term by term F = —dS / da we find the first term in optical expansion of the Casimir 
energy S (the arbitrary constant is chosen so that S —* when a oo) as 

awcos'^6' 

^iu+3u = z --2 ni.io 

27r2(4a2 -i/;W^) 

and from this one obtains the torque 

^ 2aw (w^ — 4a2) cos^ 6'sin6' 

^iu+3u = — -. — ^-3 — (ni.ii) 

7r2(4a2 — w'^sin 9) 

Analogously we can calculate the contribution to the pressure P of the two reflections paths 
2u and 2d. Again the contributions of the two paths are identical and the result simplifies 
to 

2 1 1 

and using £2 = 2^/x^~+~z^ sinO we find 

p 1 

IQn^sm^ex^ 

which integrated from Xm = (a/ sin 6* — w/2)/ cos6 and xm = (0/ sin6' + w/2)/ cos6 gives 



(III.13) 



the force along the z axis due to these paths: 



cos^^ / 1 1 



2«+2'i 487r2sin2^ V(a/sin^- w/2)3 {a/sin9 + w/2y 
This expression can now be expanded for ^ ^ 1 (quasi-parallel plates) 



(III.14) 



2u+2d 



1 



167r2 6a6 
13 



Notice that this expression vanishes when 5 — 0, as it should since for parallel plates all the 
contributions of even reflections paths vanish. 

The next term in the series is i^sd+sd, whose calculation is performed in the same fashion. 
The result is: 

3 cos5 2e /I 1 



3d+5d 



167r2 sin^ 29 \ {a/ sin 6 - w/2f (a/ sin 6 + w/2f 



(III.16) 



1 (^^^n^-^^,^ ^^^^^^^^ 



167r2 \lQa^ 32a' 
We can also present the term given by the 4 reflections paths, 



cos^ 26 f 1 1 \ I / w 



iu+id 



487r2 sin^ 29 \ (a/ sin 9 - w/2f (a/ sin 9 + w/2f J IGtt^ \Aa^ 

(III.18) 



The terms independent of 9 can be seen to reconstruct the parallel limit case F = — (1 + 
1/16 + 1/81 + ...)3/167rV. 

n 

Term by term, this series for the force reproduces the series in Ref. |lll]. The series for 
the energy and the torque agree as well. The results of the pressure method then coincide 
with those of the energy method (as for all the examples analyzed in this paper). In Ref. 
we discussed at some length the predictions of the optical method for the Casimir torsion 
pendulum. We will not repeat them here, referring the reader to that paper for further 
details. 



C. Sphere and Plane 

The sphere facing a plane is an important example for several reasons: it has been 
analyzed theoretically with various exact or approximate numerical techniques it is 

an experimentally relevant configuration; the exact solution is unknown and probably will 
escape analytical methods for a long time to come. We have already calculated the optical 
approximation to the Casimir energy in Ref. up to 5 reflections. In this paper we 
study this problem for mainly pedagogical purposes, leaving a more accurate and complete 
numerical analysis for the future. We believe it is worth studying this example because, 
contrary to the previous two examples, the enlargement factor plays an important role and 
moreover we will reanalyze this example with finite temperature in Section IV B 2. 

We calculate the pressure (and by integrating, the force) exerted on the plate by the 
sphere which, of course, equals the force exerted by the plate on the sphere. We start 
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from the qualitative observation that the rank of the contributions is the same as in the 
parallel plates case in the limit a/R 0. In all the examples we have analyzed this rank 
is preserved for any value of a/R. Moreover the ratios of the contributions to the force 
^3+5(0, -R)/Fi+3(a, i?), F4^{a, R) / F2{a, R) etc. decrease quickly as a/R increases, we believe 
due to the growing importance of the enlargement factor. 

In this paper we calculate analytically the Is term (here s stands for 'sphere' and p for 
'plate') and by using the relation Pis+ss = Pu + Pss = 2Pis proved in Section |Tll (the notation 
is the same as in that section) we are able to include the 3s term as well. 

Using the expressions for the length and enlargement factor for the Is path obtained in 
Ref. we get 



R 



[p-^ia + R-zf + p'^^ {{a + R-zf + p") ^'^ 



327r2 dz^ 



III.19) 



The final expression for the pressure Pu+^s obtained after the derivatives are taken is rather 
long, however the contribution to the force on the plate, Fu+'is (obtained by integration of 
Pis+zs over the infinite plate) is quite simple: 

/■°° hcR 
Fu+^s = 27r / dppPu+-is = (111.20) 
Jo Svra^ 

This is the largest of the contributions and increasing a/R improves the convergence of the 
series due to the presence of the enlargement factor, so the asymptotic behavior at large 
a/R predicted by the optical approximation is that given by this formula, i.e. F oc R/a^ 
or E (X R/a^. This asymptotic law is in accordance with the numerics of Ref. jl^ and the 
predictions of other semiclassical methods |18] . However, eq. ()III.20|) is in disagreement with 



the Casimir-Polder law 



2l! | which predicts E oc R? /a'^ for a^ R. This is no great surprise. 



since our method is not valid for a/R^ 1, the semiclassical reflections being corrected and 
eventually overshadowed by diffractive contributions |3, • 

We have calculated the contribution of the two reflections paths analytically as well. The 
calculation is more involved than the one reflection term but a big simplification occurs if 
one notices that, for the purpose of taking the second derivative with respect to 2; at 2; = 0, 
one can leave the reflection point on the sphere fixed. We could not prove a similar result 
for any other reflection. It is certainly not true for odd reflections but one can conjecture it 
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FIG. 2: The magnitude of the total pressure up to reflection bp in units of hc/R^ as a function 
of the radial coordinate on the plate, p/R. Upward, or red to blue a/R = 1, a/R = 0.1 and 
a/R = 0.01. 

to be true for even reflections. In this paper we have not calculated the 4 reflection terms 
and hence we could not check this conjecture for more than 2 reflections. 

And finally, we have calculated the 3p (or sps) and hence obtained the 5p, or pspsp, paths 
contribution P^p+^p] Psp+^p in the parallel plates limit should account for ~ 1/16 of the total 
force. This contribution, unlike the previous ones, must be calculated partly numerically, 
mainly because finding the refiection point on the sphere requires the (unique) solution of a 
transcendental equation. This task is achieved much more quickly by a numerical algorithm 
than by patching together the several branches of the analytic solution. 

The total pressure is plotted in Fig. in |21 while the various contributions (keeping in mind 
that Pis+3s and Psp+^p are negative and P2+2 is mainly positive) are shown in Fig. El Figure|21 
reveals some interesting features of the pressure in this geometry: the total pressure decays 
very quickly with the distance as P ~ the exponent a seems to depend upon the 

distance a/R, but for a/R < 0.1 a good fit is obtained with a = 6, in accordance with the 
asymptotic expansion of the 1 + 3 refiection term Eq. piI.19|l : by decreasing the distance 
between the sphere and the plate, the pressure becomes more and more concentrated near 
the tip, giving us reasons to trust our approximation and supporting the use of the PFA 
as a first approximation in the limit a/R 0. Figure El shows the relative importance 
of the contributions due to the different paths. As expected the contribution to the total 
pressure decreases quite fast by increasing the number of refiections. In Fig. |31 one can also 
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FIG. 3: Contributions to the pressure in units of hc/R^ as a function of p/R, for fixed a/R = 0.1. 
Downward or red to blue, we have — Pi^+ss, — Psp+sp and P2+2- Although unnoticeable in this 
figure, the curve P2+2 changes sign at around p/R ~ 0.4 (see Figure H] for a similar situation). 
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FIG. 4: Contribution of the two reflection path(s) to the pressure in units of hc/R'^ as a function 
of p/R, for fixed a/R = 0.01. The pressure becomes negative, showing that the sign of the pressure 
is not determined by the number of reflection only. 

see that the sign of the pressure is not determined simply by the number of reflections of 
the underlying optical path — as for the contribution to the energy density. 

By integrating the pressure over the whole plate we obtain the force F. It is useful to 
factor out the most divergent term of the force, as predicted by the PFA, so we deflne the 
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FIG. 5: The ratio between the optical force up to the 5p reflection and the most divergent term 
in the PFA, as defined by eq. pil.21|) . 



quantity f{a/R) as 



F{a) 



TT^R 



f{a/R). 



(III.21) 



720a3' 

Since we include only a finite number of reflections it is convenient to factor out the constant 
C(4)/(l + 1/16) such that / is normalized with /(O) = 1. The function f{a/R), calculated 
including paths Is, 3s, 2, 3p and 5p, is plotted in Figure El When a/ R f is fitted by 



f{a/R) = 1 - 0.10 a/R + O {{a/Rf) . 
By comparing to the results of [l^ 

hner^{a/R) = 1 + 0.05 a/R + O {{a/Rf) 



(III.22) 



(III.23) 



there is the difference in the sub-leading term. 

By neglecting the 5s + 7p reflection paths (which in the parallel plates case contribute 
~ 2% of the total force) we can only assert that the functions / in pil.22|l and pil.23jl 
represent the optical approximation with an error of 2%. When plotted on the whole range 
of a/R where the optical approximation is to be trusted the pressure and energy method 
curves never differ more than 2%. However there is no such a bound on the sub-leading 
term which, on the contary, depends on the higher reflections contributions which have not 
been included in this calculations.^ With the terms calculated at this point, we cannot 



For example consider that including only Is, 3p, and 2 and reflections would have given a sub-leading 
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make a precise statement about the sub-leading term. We can however safely say that 
the subleading term a/R coefficient is quite small and our method disagrees with the PFA 
prediction —0.5a/R. The sphere opposite plate is such an experimentally relevant geometry 
that further, more accurate studies need to be performed to compare with experimental 
data. 

In conclusion, the lessons to be learned from this example are two: 1) The calculations 
with the pressure method are even quicker and simpler than the energy method and some- 
times can give analytic results for non-trivial geometries and 2) the sub-leading terms must 
be compared only between calculations performed with the same accuracy.^ 



IV. CASIMIR THERMODYNAMICS. 



As measurements of Casimir forces increase in accuracy they become sensitive to thermal 
effects. The natural scale for Casimir thermodynamics is a distance, /3 = hc/irT, which at 
room temperature is about 2.5 microns. [To avoid confusion with the wave number k, we 
set Boltzmann's constant equal to unity and measure temperature in units of energy. We 
continue to keep h and c explicit.] So, assuming the corrections are of O ^(a/^)"j, depend- 
ing on the value of a thermal effects might be expected between the 10% (for a = 1) and 
0.3% (for a = 4, the standard parallel plates result) level for Casimir force measurements on 
the micron scale. In open geometries, like the sphere and plane, even longer distance scales 
are probed by Casimir effects, and this gives rise to interesting changes in the temperature 
dependence of the Casimir free energy in comparison with the case of parallel plates P|. The 
optical approximation is well suited for discussion of thermodynamics since the thermody- 
namic observables, like the Casimir energy, can be expressed in terms of the propagator. 
Here we consider again a non-interacting, scalar field outside rigid bodies on which it obeys 
Dirichlet boundary conditions. 

term —0.16a/R instead of —O.lOa/R in Eq. (|III.22|I . The sub-leading term then changes of 50% by adding 
the 3s + 5p reflection terms which contributes only up to 8% of the total. 

AS would like to thank M. Schaden and S. Fulling for conversations on this point during the workshop 
'Semiclassical Approximations to Vacuum Energy' held at Texas A & M, College Station, TX, January 
2005. The concerns about the errors to be associated with the optical, semiclassical or proximity force 
approximation is still open to debate and is strictly connected to one of the most challenging open 
problems in spectral theory i. e. how to go beyond the semiclassical approximation to the density of states 
of a positive Hermitian operator. 
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Before entering into a technical discussion of temperature effects, it is useful to anticipate 
one of our central results which follows from qualitative observations alone. As T ^ the 
temperature effects probe ever longer distances. Even at room temperature the natural 
thermal scale is an order of magnitude larger than the separation between the surfaces in 
present experiments (see Ref. |^]). Since long paths contribute little to the Casimir force, 
we can be confident that thermal effects vanish quickly at low temperature. However, the 
leading T-dependence at small T comes from regions beyond the range of validity of the 
optical (or any other) approximation, so we are unable to say definitively how they vanish 
for geometries where no exact solution is possible {i.e. other than infinite parallel plates). 

This section is organized as follows: First we discuss the free energy and check our 
methods on the parallel plates geometry; then we discuss the temperature dependence of 
the pressure, which we apply to the sphere and plate case. Finally we discuss the difficulties 
associated with the T limit. 

A. Free Energy 

The free energy is all one needs to calculate both thermodynamic corrections to the 
Casimir force and Casimir contributions to thermodynamic properties like the specific heat 
and pressure. However like the Casimir energy, Casimir contributions to the specific heat, 
pressure, etc., are cutoff dependent and cannot be defined (or measured) independent of 
the materials which make up the full system. So we confine ourselves here to the thermal 
corrections to the Casimir force. The problem of parallel plates has been addressed before 
and our results agree with those j^. 

1. Derivation 

We start from the expression of the free energy for the scalar field as a sum over modes 




= 5^ In (1 - e-'3('^"-'^)) + J2 



n n 



= J^ + S, 



(IV.l) 
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where n is the chemical potential, and the last term is the Casimir energy, or the free energy 
at zero temperature, since JF = for T = 0. The Casimir energy S, being independent of 
the temperature, does not contribute to the thermodynamic properties of the system. It 
however does contribute to the pressures and forces between two bodies. The force between 
two bodies, say a and b, is obtained by taking the gradient of the free energy with respect 
to the relative distance fab 

fab = -Vab:F. (IV.2) 

At T = we recover the familiar result / = — V£^. 

Next we turn the sum over modes into a sum over optical paths. Following the same 
steps that led from Eq. pi.4|) to Eq. ()11.15|) we obtain 

//»oo 
d^x J dk p{x,k)\n{l-e-^^'^^^^-^^) . (IV.3) 

where p{x, k) is given by Eq. pi.l6j) . By specializing to a massless field in 3 dimensions with 
zero chemical potential (to mimic the photon field), and substituting the optical approxi- 
mation for the propagator Eq. ()II.13|) . we obtain the sum over paths 

^ = V = y2{-lY-\- [ d^x r dk k sm{kir) In (l - e"^'^^) . (IV.4) 

r=0 r ' •^'^r Jo 

Here the term JFq, the direct path, gives the usual free energy for scalar black body radiation. 
Using the values for the direct path, we have Aq = I/^q ^"^^ = — a;| ^0 when taking 
x' — >■ X. We get the familiar textbook expression 

where V is the (possibly infinite) volume outside the bodies. 

The general term JF^ associated with the path r is calculated by performing the k integral 
in Eq. (jEl): 

= f rf3^^i/2 1 r_2 + I. [coth^ + ^csch^f,) 1 (IV.6) 

where ir = irT^T/hc = £r/(3 measures the path length relative to the thermal length scale. 

Eq. ()IV.6|) is the fundamental result of this section and gives a simple, approximate 
description of thermal Casimir effects for geometries where diffraction is not too important. 
There are no divergences in any of the J?>, ultraviolet or otherwise, even for the direct path 
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(as we saw in eq. ()IV.5|) ) and the first reflection path. All the ultraviolet divergences are 
contained in the Casimir energy £. Indeed, by expanding the integrand of equation pV.6|) 
at short distances, i.e. -C 1, we obtain 



' m L 



-2 + 4 coth4 + 4csch^ 



/?3 



(73 _|_ 



(IV.7) 



.45/? 945/?3 

Only the 1-reflection path length can go to zero to generate a divergence. For this contribu- 
tion Ar diverges like as — ^ 0, however this is compensated by the ir term in pV.7|l 
so the expression is finite and then integrable. 

To check for infrared divergences notice that at large distances, £r ^ 1, the integrand of 
PV.6|) goes to ~ A^^/f^. For an infinite flat plate the A^ ~ 1/^;^, where z is the normal 
coordinate to the plate, and the integral is hence ~ dz/z^ at large z. For finite plates the 
domain of integration is finite and for curved plates the enlargement factor falls even faster 
than 1/^^, and the integral remains convergent. 

Since the integral converges in both the infrared and ultraviolet, it is safe to estimate the 
important regions of integration by naive dimensional analysis. This leads to the conclusion 
that The paths that dominate the temperature dependence of the Casimir force have lengths of 
order the thermal length (3. High temperature implies short paths. Very low temperatures are 
sensitive to very long paths. Long paths involve both paths experiencing many reflections, 
which are sensitive to the actual dynamics at and inside the metallic surface, or paths 
making long excursions in an open geometry, which are sensitive to diffraction. Either way, 
low temperatures will present a challenge. 



2. Parallel Plates 



We know that in the limit of infinite, parallel plates the optical approximation to the 
propagator becomes exact. Hence our method gives another way to calculate the free energy 
of this configuration of conductors. It is convenient to study this example to check against 
known results and to prepare the way for a study of the T — ^> limit. 

We recall that for this configuration the expression for the enlargement factor is A = 1/ ^^ 
and the lengths are given by = "^na (where a is the distance between the plates) and 
^2n+i, u = 2(a — 2;) + 2na, i2n+i, d = 2^ + 2na, the notation being the same as in Section 
IIII Al should at this point be familiar to the reader. 
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As in the zero temperature case it is useful consider even and odd reflection contributions 
separately and as for the zero temperature case, the sum over odd reflections turns into an 
integral over z from to oo 

-^odd = y2 J^2n+i,d + J^2n+i,u = — ^^•S' / dx— [-2 + a;(cothx + X csch\)] , (IV.8) 

n=0 '^^ P 

where x = 2z/ (3 and S is the area of the plate. The definite integral can be easily performed 
numerically and its value is u = 0.06089..., 

which is independent of the separation, a, and therefore does not contribute to the force. 

Let us turn now to the even reflection paths. They have constant length 2na, so the 
volume integral simply yields the volume between the surfaces v = Sa. We already calcu- 
lated the zero-reflection term jFg in Eq. pV.5|) . The remaining even reflection contributions 
(2,4,6,... reflections) J^even,r>2 can be written as an infinite sum 

he 1 °° 1 

^even,r>2 = -2 — Sa^ ^ — [-2 + X„ (cothx„ + X„CSch^X„)] (IV.IO) 

P n=l 

where x„ = 2na/ (3 = nr (this defines the dimensionless temperature r) and we have intro- 
duced an overall factor of two to take into account the multiplicity of the paths. Thus the 
total free energy for parallel plates is the sum of jFg (eq. ()IV.5|) and the results of eqs. pV.9|) 
and (HVlOll ). 

= + ^,Su - -^Sa E ^ + (cotha;„ + x.csch^x^)] (IV.ll) 

It is not possible to rewrite jFy in a closed form, but the sum is easy to compute numerically 
and the high and low temperature expansions are easy to obtain analytically. At high 
temperatures (and fixed a) r ^ cxd, and the summand gin) in eq. pV.ll|l falls rapidly 
enough with n 

9{n) = [-2 + (rn) (coth(™) + (™)csch2(™))] = -1- [-2 + rn] + O (e"™) , 

(IV.12) 

that the limit may be taken under the summation, with the result. 



^ 3ven,r>2 — n '^'^ / . 

" n=l 



CO r ^ 1 



2^3^3 ^4^4 
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'^'\ST,^,S. (IV,13) 
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Notice that the second term cancels the even paths contribution to the Casimir energy. 
Hence the final expression for the high T expansion of the free energy is particularly simple, 

The first term is usual black body contribution to the bulk free energy. It does not contribute 
to the force. The second term is also independent of a and does not give rise to any force. 
The third term instead gives the thermal Casimir force. Notice that he has disappeared 
from this expression. Called the "classical limit" , this high temperature behavior has been 
noted before and some early results are even due to Einstein (in pg. 2; see also In 
the next section, after the thermal corrections to the pressure are calculated, we show how 
to extend this result to other geometries. 

Note some interesting features of the T — ^ oo limit: First, the sum over paths converges 
like the sum of as indicated by the appearance of C(3)- While slower than the T = 

convergence, it is still rapid enough to obtain a good approximation from low refiections. 
Second, note that the T —>■ oo problem in 3-dimensions corresponds exactly to a T = 
problem in 2-dimensions. This is an example of the familiar dimensional reduction expected 
as T ^ oo. We can give a short proof of this result. Let us first write: 

F = -hogZ (IV.15) 

P 

where Z is the partition function. We need to evaluate Z to the lowest order in (3 when 
j3 0. The thermal scalar field theory can be written as a free theory on the cylinder 
M'^ X [0,(3). For /5 — the dynamics along the thermal coordinate is frozen in the ground 
state, with energy Eq = 0, where does not depend on the thermal coordinate. The partition 
function Z is now Z = Z3 + O (e^^^^) where Ei is the first excited state Ei oc and Z3 is 
the partition function of the remaining three-dimensional problem in M^. If the conductors 
geometry is symmetric along one spatial coordinate, say x (in the parallel plates problem 
we have two of these directions, x and y) this can now be interpreted as an Euclideanized 
time variable extending from to L^/c. So we will write Z^ = Z2+1 = e~^^^^'^^'^ where £2 is 
the Casimir energy of the 2 dimensional problem of two lines of length Ly, distant a. The 
free energy F is then: 



~\ogZ ^ -ilogZ2+i = T^^S2 = -TL^Ly^ 
jj jj he IGn^a 



24 



Since 5" = L^Ly This is exactly the a-dependent term in eq. f)IV.14|) . If the geometry is not 
translational invariant then we can only say from eq. pV.16|) that the free energy is linear 
in T (since Z2+1 is independent of /?). Later, by using the optical approximation we will 
find an explicit analytic expression valid also for non- symmetric, smooth geometries. 

For low temperatures, r — 0, the terms in the n-sum in eq. fjlV.ll|) differ very little from 



each other so we can use the Euler-McLaurin formula 



23, 



/"OO 11 7/1 

J1^9{n) = dx g{x) - -g{Q) - -g\Q) + ... =--- + 0{t). (IV.17) 

Substituting into eq. pV.lljl we find that the first term in eq. pV.17jl cancels the sum over 
odd reflections (the second term in eq. ()1V.11|) ) and that the second term in eq. ()1V.17|) 
combines with JFq to give a very simple result, 

- ^ " 90(^3 ■ (IV.18) 

at low temperatures. This has a simple physical interpretation: the typical thermal excita- 
tions of the field at low temperature have very long wavelengths, it is hence energetically 
inconvenient for them to live between the two plates. As a result the only modification of the 
T = result is to exclude from the standard black body free energy the contribution from 
the volume between the plates. One could imagine measuring this effect as a diminished 
heat capacity for a stack of conducting plates inside a cavity. 

The low temperature result, eq. pV.18|) . is deceptively simple. Its simplicity obscures an 
underlying problem with the T — > limit. We postpone further discussion until we have 
explored the temperature dependence of the pressure. Suffice it to say for the moment, that 
eq. pV.lSj) probably does not apply to realistic conductor with finite absorption, surface 
roughness, and other non-ideal characteristics. 



B. Temperature dependence of the pressure 

In this section we will obtain the temperature dependence of the pressure within our 
approximation and apply it to a preliminary study of the sphere and plate case. To begin, 
we calculate the thermal average of an operator O quadratic in the real scalar field 0. The 
average of a generic operator O is given by the trace over a complete set of eigenstates \^q) 
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of the Hamiltonian weighted by a Boltzmann factor: 



{O)^ = 5^e-^^'^(^JO|^„). (IV.19) 

a 

After some algebra we find 

j 

j 

where denotes the thermal average, j labels the normal modes ipj (cf. Section II), nj is 
the occupation number of the mode j and Ej its energy. The quantities Oj are read from the 
decomposition of the diagonal part of the operator O written as Cdiag = Cj(a]aj + ajaj) 
where aj is the annihilation operator of the mode j. 

The Oj for the pressure can be read easily from the analysis in Section |n] 

P,= Mm ^^L^nMx')^P,{x) (IV.21) 

So we can write the pressure on the plate at non-zero temperature as 



f°° k /I 



where we have used Eq. pi.30|) . 

Next we introduce the optical approximation for the propagator and limit ourselves to 
massless scalars E{k) = hck. The discussion of the divergences parallels that of Section HTl 
and needs not be repeated here. We remove Pq and Pi and leave all the finite contributions 
r eTZ. The optical approximation for the pressure exerted by a massless scalar field reads 



A 



g[icoth(^,/,§) 



(IV.23) 



where it is understood that the zeroth and first reflection terms, which contribute to the 
pressure on each surface individually, but not to the force between surfaces, have been 
dropped. 
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Before applying this to the sphere and plate problem, let us again look at the limiting 
behavior as T ^ oo and T — » 0, and draw some conclusions independent of the detailed 
geometry. First consider T — oo. The shortest paths in the sum in eq. ()1V.23|) are of order a, 
the intersurface separation. [Remember that the optical approximation is accurate as long as 
the important paths are short compared to R, a typical radius of curvature of the surfaces.] 
At high T we can take the /? — limit under the sum over reflections since the resulting sum 
still converges. Therefore low reflections dominate, and we can see, retrospectively, that the 
high temperature approximation applies when /3/a— >0. SoasT— >oo, 

1/2 

(IV.24) 



This limit has been called (it has been previously found for the parallel plates case) the 
"classical limit" 0, Q, [2^ , since the final expression for high temperatures, reinserting h 
and c, 

a1/2 

P ^ (IV-25) 

r 

is independent of h and c apart from exponentially small terms. This expression amounts in 
neglecting the 1 in the expression {2nj + 1)^, corresponding to normal ordering or neglecting 
the contribution of the vacuum state. 

At low temperatures, /? ^ 00, it is not possible to interchange the limit with the sum. 
The relevant quantity is i coth(£r//3), which goes like 

icoth (^)=^ + -k ^1^ + (S-) (IV.26) 

P \/3j ir 3/32 45/34 

as /3 00. The first term yields the familiar T = expression. The others would give 
divergent contributions because of the factors of ir in the numerators (even after the in- 
clusion of the enlargement factor A^). Of course the sum over reflections of the difference, 
i coth(£r//?) — J-, converges to zero as /? — 00, so thermal corrections definitely vanish for 
any geometry as T — as expected. Once again we relegate more detailed consideration of 
the T ^ limit to a later subsection. 



1. Sphere and plate 



In this section we calculate the pressure and total force for the configuration of a sphere 
facing a plane at non-zero temperature within 5p reflections. The optical approximation 
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FIG. 6: The p dependence of the Is + 3s contribution to the pressure Pis+zs for the sphere and 
the plate in units of hc/R^ for various temperatures. Two effects must be noticed. The top 3 
curves (in blue) show the high-temperature region where the pressure is proportional to T (notice 
the logarithmic scale). The two lower curves (in orange and red) show the low-temperature region 
when increasing the temperature changes the asymptotic behavior of P for large p {i.e. p ^ 13) 
while for small p the behavior reduces to the zero-temperature limit. 

should be accurate if the important paths are short compared to R, the radius of the sphere. 
On the other hand the thermal corrections to the force are sensitive to paths with lengths 
of order (3. So we must have R ^ (3 and i? ^ a in order to obtain reliable results from 
the optical approximation. Fortunately this is a region of experimental interest: present 
experiments use, for example, a ^ 0.5fim, R ^ lOO/im, and at room temperature, f3 ~ 
2.5/im. In this regime the optical approximation should give a good description of the 
thermal corrections to the force between perfectly reflective, perfectly smooth conductors. 

The expression for the pressure is given by Eq. pV.23j) . the enlargement factors and 
lengths are the same as in the T = case. By applying Eq. ()1V.23|) to the Is + 3s paths we 
find the results in Figure IHl Notice that at high temperatures increasing the temperature 
essentially scales the whole plot proportionally to T. The force is then linearly dependent 
on the temperature (this is the 'classical limit' already discussed in Section IIVB|) . More 
details are given in the caption of Figure El 

A dimensionless function f{a/R,P/R) can again be defined by rescaling the total force 
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F to extract the leading term as a — > 0. The hmiting behavior a ^ is not affected by 
temperature effects so we stick to the old definition for /: 

In Figure d we present / (up to 5 reflections) for 5 different values of P/R (we choose 1, 
1/2, 1/4, 1/8 and 1/16 recognizing that /5 ~ 1 strains the limits of our approximations) 
and varying a. Notice that in a neighborhood oi a/R = 0, shrinking as (3/R increases, the 
function / is very well approximated by the T = form, already discussed in Section Fill C\ 
f{a/R) ^ 1 — O.la/R. It is not useful to study the derivative A{P/R) = df{x,P)/dx as 
X = a/R ^ since this will take the constant value predicted by the zero temperature 
analysis, or —0.1 in this approximation, for any value of the temperature we choose. 

It is also clear from the previous discussions leading to equation pV.24|) that in the 
opposite regime, for a//3 1, we must have F oc R/a?(3 = RT/a? (the 'classical limit'). In 
fact, the first term in the high temperature expansion pV.24|l integrated over p converges 
and gives a finite force linear in T. For this problem, the first term in the reflection expansion 
for high temperatures can even be calculated analytically: 

R . R 

Unfortunately there is no such simple closed expression for higher reflection terms (nor for 
this first term at arbitrary T). However, if one believes that the rank of contributions is 
similar to the parallel plates case one should feel safe to say that this truncation captures 
the optical approximation within a C(3) — 1 — 20%. Hence our statements are at least 
qualitatively correct. 

This expression for the force gives a prediction for the function /, defined in Eq. pV.27|) . 
At this level of accuracy (Is + 3s reflection) and for a//3 > 1, apart for exponentially small 
terms in the temperature expansion we have 

90 a i? . X 

/,.«. - (IV.29) 

which grows linearly in a/R and is (interestingly enough) independent of R. This is evident 
in Fig. Elfor the curves with j3 = 1/8, 1/16. For higher (3 the linear growth starts at higher 
values of a not shown in Fig. [7| Moreover the exponential accuracy manifests itself in the 
sudden change of behavior from / ~ 1 — O.la/R to / cx a/p. 
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Fu+ss = -hc-^ + O (e-^/^) ^ -nc At. (IV.28) 



0.02 0.04 0.06 0.08 0.1 0.12 
a/R 



FIG. 7: The function f{a/R, 13/R) as a function of a/R for P/R (from red to violet or down 
up) = 1, 1/2, 1/4, 1/8, 1/16. /(O) ~ 0.98 since we are summing only up to reflection 5p. The two 
lowermost curves, red and orange (/? = 1, 1/2) superpose almost exactly. 

It is quite easy to extract a universal prediction from this data, whatever the definitive 
numbers are, after the sum over optical paths is carried to sufficiently high order: for any 
non-zero temperature the function f{a/R) will deviate from his zero-temperature behavior 
at a ^ j3 ^ hc/T. The deviation will be in the upward direction, increasing the attractive 
force between the bodies. Eventually, for sufficiently large distances, the high temperature 
behavior given by eq. pV.25|) (or pV.29|) for the sphere-plane problem) will be recovered. 

C. Thermal corrections at low temperatures 

The preceding examples have made it clear that in the language of the optical approx- 
imation, thermal corrections at low temperature arise from very long paths, ir ~ This 
can be seen from the general form of the free energy, eq. pV.6|) . or in the attempt to take 
the (3 oo limit under the summation in eq. pV.23j) . which fails because of the expansion, 
eq. pV.26|l . Here we examine this non- uniformity more carefully in general and in particular 
for the parallel plate case, where all the expressions are available. We then attempt to draw 
some conclusions about the magnitude of corrections at low temperature and the possibility 
of calculating them reliably in an model that idealizes the behavior of materials. 

We return to eq. ()IV.22|) . which gives the exact expression for the pressure, and separate 
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out the thermal contribution, 




(IV.30) 



still exact. Expanding the denominator in a geometric series, we find 




—mfihck 



(IV.31) 



Each term in the sum is a Laplace transform of the Greens function. Clearly, as /? ^ oo the 
frequencies that dominate this integral are oc ~ T. 

What are the low frequency contributions to Q{x',x,k)l In the ideal case of infinite. 



However in realistic situations the plates arc finite and/or curved, the geometry is open, 
and there is no gap in the spectrum. The low-Zc part of the spectrum is sensitive to the 
global geometry, including edges and curvature, and to the low frequency properties of the 
material. If the conditions are close to the ideal, the contributions to 5P from small k may 
be small. However as T — > 0, they dominate. We conclude that the T — > behavior of 6P 
cannot be calculated for realistic situations. 

The optical approximation does not take account of diff'raction, and cannot accurately 
describe the T — > limit. Nevertheless it is interesting to see how it fails, since this sheds 
light on the problem in general. Substituting the optical expansion for the Greens function 
(replacing 9^,^ ~^ ^'^^ setting h — c—1) v^e find 



The problems with T — > are quite apparent: as /3 — > oo all paths become important. 

Next we specialize to parallel plates where Ir — (2ar±22;). The derivative can be carried 
out explicitly. For simplicity we focus on m = 1 {5P — X]m=i ^Pm)-, 



perfectly conducting, parallel plates, there is a gap in the spectrum at low k: k > ^. 





(IV.32) 




(IV.33) 



which can be rewritten using the variable r = 27ra//3 introduced earlier. 




(IV.34) 
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The sum can be performed, giving 

1 



5Pi = - 



1 



coth 



71(3 

2a 



csch 



71/3 

2a 



(IV.35) 



The second term in brackets is exponentially small as (3 
m-dependence, and sum over m, we obtain 

^2 



6P 



oo. If we ignore it, restore the 



(IV.36) 



which agrees with our earlier calculation, as it must. 

However eq. pV.34j) allows us to study the convergence of the sum over reflections as 
(3 — > oo. Instead of performing the sum analytically, we sum up to some rmax = X- Since 
r ^ 0, we can once again use Euler-Maclaurin, to rewrite the sum over r as 



5Pi = - 



2 1 
^]3^ 



X 



+ 



1 3r2xV7r2- 1 



(IV.37) 



_2 (l + r2XV7r2)2 2 (1 + r2XV7r2)3 ' 

where the omitted terms are higher Euler-Maclaurin contributions that are unimportant as 
P ^ oo {i.e. T 0). 

If the upper limit on the sum, X, is taken to oo, only the first term, 1/2, survives and 
gives the expected result. The question is: How large must X be before the limiting behavior 
set in? Dropping the third term in eq. ()IV.37|) . which is subdominant, we can rewrite 6Pi 
as 



X 



;i + X2r2/7r2)2 



7r2/?4 



-firX) 



(IV.38) 



The function f{z) is negative definite and has a minimum at z = 1/2^3 ^ 0.29 where 
it takes the value —3^/^/32 ~ —0.16. So in order the result Eq. pV.36|) to be valid we 
must include X ^ X^ = ti/VSt terms in the sum. For example in a typical experimental 
situation we have a = 0.5 fim and T = 3007^^ so /5 = 8/im, r = S/vr and X^ = S/Vs = 4.6. 
In this case it is necessary to go to X ~ 20 before the contribution of |/(rX)/r| is smaller 
than 1/2. This means paths with ~ 40 reflections and path lengths of order 20 fim. With 
40 chances to sample the surface dynamics of the material and paths of 20/im available to 
wander away from the parallel plate regime, the idealizations behind the standard parallel 
plates calculation must be called into question. 

It must be said however that in the modern experiments the temperature corrections are 
at most of the order of a few percent at a ~ Ifim and vanish when a —>■ 0. Nonetheless 
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we want to point out that there is a conceptual difference between formulations based on 
the infinite parallel plates approximation, extended to curved geometries by means of the 
PFA, and a derivation (like ours) in which the curvature is inserted ab initio. The thermal 

nn 

and curvature scales interplay in a way that the usual derivations [3|, y| could not possibly 
capture, giving rise to different power law corrections in a/ (3. It is worth reminding the 
reader that the usual numerical estimates of thermal corrections are based on the infinite 
parallel plates power law (a//3)^. A smaller power like {a/ (3)"^ would give a much bigger 
upper bound. 

To summarize: temperature corrections are small at small T, but the existing methods 
of calculating them, including both our optical approximation and the traditional parallel 
plates idealization, cannot be trusted to give a reliable estimate of the T-dependence at 
small T. 



V. CONCLUSIONS 



In this paper we have shown how to adapt the optical approximation to the study of 
local observables. We have illustrated the method by studying the pressure, but the method 
applies as well to other components of the stress tensor, to charge densities, or any quantity 
that can be written in terms of the single particle Greens function. The advantage of the 
optical approximation is to extend the study of these local observables to novel geometries. 
In particular we developed an expression for the Casimir pressure on the bodies and applied 
our main result Eq. pi.31|) to the study of three important examples: parallel plates, the 
Casimir pendulum and a sphere opposite a plate. 

We have also shown how to calculate, within this approximation scheme, thermodynamic 
quantities and thermal corrections to the pressure in the general case and applied our results 
to the example of parallel plates (retrieving the known results) and to the case of a sphere 
opposite a plate. Along the way we have given a proof of the "classical limit" of Casimir 
force for any geometry (within our approximation), i.e. the fact that Casimir forces at 
high temperatures are proportional to the temperature and independent of h, a fact that 
previously was known only for parallel plates. 

Finally, we argued that all known methods of computing the temperature dependence of 
the Casimir effect are suspect as T — > 0. 
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